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Fourier  theory  of  thermal  transport  considers  heat  transport 
as  a  diffusive  process  where  energy  flow  is  driven  by  a 
temperature  gradient.  However,  this  is  not  valid  at  length 
scales  smaller  than  the  mean  free  path  for  the  energy  carriers  in 
a  material,  which  can  be  hundreds  of  nanometres  in  crystalline 
materials  at  room  temperature.  In  this  case,  heat  flow  will 
become  'ballistic'— driven  by  direct  point-to-point  transport 
of  energy  quanta1.  Past  experiments  have  demonstrated  size- 
dependent  ballistic  thermal  transport  through  nanostructures 
such  as  thin  films,  superlattices,  nanowires  and  carbon 
nanotubes1-8.  The  Fourier  law  should  also  break  down  in  the 
case  of  heat  dissipation  from  a  nanoscale  heat  source  into  the 
bulk.  However,  despite  considerable  theoretical  discussion  and 
direct  application  to  thermal  management  in  nanoelectronics2, 
nano-enabled  energy  systems9,10  and  nanomedicine11,  this 
non-Fourier  heat  dissipation  has  not  been  experimentally 
observed  so  far.  Here,  we  report  the  first  observation  and 
quantitative  measurements  of  this  transition  from  diffusive  to 
ballistic  thermal  transport  from  a  nanoscale  hotspot,  finding 
a  significant  (as  much  as  three  times)  decrease  in  energy 
transport  away  from  the  nanoscale  heat  source  compared  with 
Fourier-law  predictions. 

When  heat  flows  from  a  hotspot  with  a  dimension  L  smaller  than 
the  mean  free  path  of  phonons  in  the  heat  sink  A ,  energy- carrying 
phonons  travel  ballistically  away  from  the  source  for  a  significant 
distance  before  experiencing  collisions.  However,  in  a  Fourier 
picture  of  thermal  transport,  heat  flow  is  assumed  to  be  diffusive 
(experiencing  many  collisions),  with  the  thermal  flux  q  determined 
by  the  temperature  gradient  (g  ~  VT).  When  there  is  a  significant 
energy  density  change  over  a  distance  less  than  A ,  there  is  no  local 
phonon  equilibrium  and  therefore  temperature  cannot  be  defined, 
as  illustrated  in  Fig.  la.  Under  these  circumstances,  the  Fourier 
law,  which  assumes  local  thermal  equilibrium,  overestimates  the 
thermal  transport1,12,13. 

In  the  intermediate  case  when  the  hotspot  dimension  L 
is  of  the  order  of  A ,  transport  is  neither  fully  ballistic  nor 
diffusive.  The  full  Boltzmann  transport  equation  for  phonons 
must  then  be  solved  to  model  thermal  transport1 .  Unfortunately, 
an  analytical  solution  of  the  Boltzmann  transport  equation  for 
real  systems  is  usually  impractical  without  some  simplifying 
assumptions,  and  numerical  simulations  vary  widely  in  their 
prediction  of  the  strength  of  this  nanoscale  hotspot  effect14-16. 
Several  phenomenological  models  have  been  proposed  to  explain 


the  reduction  in  thermal  flux  for  a  nanoscale  heat  source,  but 
there  is  significant  variation  in  their  predictions17-19.  One  model, 
the  ballistic-diffusive  equation  model17,  provides  clear  insight  by 
separating  the  size- dependent  heat  flux  into  contributions  from 
ballistic  and  diffusive  components,  similar  to  what  is  done  in  the 
case  of  rarefied  gas20  or  electron  flow21 .  Figure  2  illustrates  the 
predictions  of  this  model  qualitatively:  at  large  length  scales,  the 
thermal  transport  is  dominated  by  the  Fourier  component,  but  as 
the  dimension  of  the  hotspot  is  reduced,  ballistic  transport  becomes 
increasingly  important.  Other  models  suggest  that  the  transport  is 
determined  by  different  effects,  such  as  non-propagating  phonon 
modes22  and  interface  effects19.  This  extensive  modelling  effort 
has,  however,  taken  place  in  an  absence  of  real  experimental 
data  with  which  to  benchmark  the  models.  Given  the  variations 
between  these  models,  quantitative  experiments  are  very  important, 
both  for  a  fundamental  understanding  of  the  basic  and  important 
phenomenon  of  heat  flow,  and  to  benchmark  models  of  nanoscale 
heat  transfer  for  nanosystem  design. 

In  many  practical  nanosystems,  an  interface  separates  a  small 
heated  region  from  a  heat-sink  substrate,  and  the  thermal  resistance 
of  this  interface  must  also  be  taken  into  account.  The  heat 
flow  from  the  heated  nanostructure,  as  shown  in  Fig.  la,  can 
then  be  characterized  by  an  effective  resistance  acting  at  the 
interface — with  contributions  from  both  a  ballistic  correction  to 
the  Fourier  law  and  the  conventional  thermal  boundary  resistance 
(TBR)  characterizing  interfacial  transport  at  bulk  length  scales. 
TBR  is  measured  in  thin  films  using  time-resolved  optical  pump- 
probe  methods  such  as  transient  thermoreflectance4,23  (TTR). 
However,  the  diffraction  limit  for  the  visible  probe  light  limits  this 
technique  to  structures  with  linewidths  greater  than  ~  1  pm — too 
large  in  scale  to  observe  ballistic  size  effects  from  a  nanoscale 
heat  source  at  room  temperature.  Using  a  specialized  geometry 
that  allowed  confined  heating  without  a  heterogeneous  material 
interface,  Sverdrup  et  al.  used  electrical  resistance  thermometry 
to  measure  steady-state  temperature  distributions  in  a  thin  5  pm 
silicon  membrane  near  a  300-nm-wide  highly  doped  silicon 
resistor.  At  temperatures  below  200  K,  they  observed  a  slight 
deviation  from  the  prediction  of  Fourier  heat  conduction22. 
However,  this  experiment  did  not  allow  for  any  systematic  and 
quantitative  comparison  between  experiment  and  theory  owing  to 
the  limitations  posed  by  the  sample  geometry.  More  recent  work 
did  not  observe  any  ballistic  size  effect  at  room  temperature24. 
It  is  worth  noting  however,  that  all  of  these  measurements  were 
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Figure  1  |  Experimental  set-up  for  measuring  thermal  transport  across  a  nanoscale  interface,  a,  Schematic  illustrating  the  difference  between  diffusive 
and  quasi-ballistic  thermal  transport  across  an  interface.  When  the  size  of  the  interface  is  less  than  the  phonon  mean  path  in  the  substrate,  the  Fourier 
diffusive  heat  equation  breaks  down,  b,  Sample  geometry  showing  infrared  laser  illumination  and  soft  X-ray  detection  scheme,  c,  Experimental 
geometry:  an  infrared  laser  beam  at  800  nm  heats  the  nanostructure,  and  a  soft  X-ray  beam  at  29  nm  monitors  the  heat  flow  from  the  nanostructure  into 
the  substrate. 


done  using  steady-state  techniques  poorly  suited  to  capturing 
short-timescale  ballistic  dynamics. 

In  this  work,  we  use  a  beam  of  ultrafast  coherent  soft 
X-rays,  at  a  wavelength  of  29  nm,  to  directly  observe  the 
cooling  dynamics  of  a  nanoscale  heat  source  into  its  bulk 
surroundings.  By  interferometrically  monitoring  displacement  in 
a  heated  nanostructure  using  diffraction  of  soft  X-ray  light,  we 
can  detect  dynamic  temperature  changes  more  sensitively  than 
would  be  possible  using  either  ultrashort  optical  pulse  probing 
or  steady-state  measurements.  This  ultrafast  short-wavelength 
transient- diffraction  methodology  represents  a  newly  developed 
experimental  capability  that  allows  us  to  map  out  nanoscale  thermal 
transport  over  a  range  of  nanostructure  sizes. 

Figure  la  illustrates  the  difference  between  diffusive  and  ballistic 
heat  transfer,  and  Fig.  lb,c  shows  the  sample  geometry  and  overall 
experimental  set-up.  Periodic  nickel  lines  with  a  height  of  20  nm, 
length  of  120  pm  and  with  widths  L  varying  from  65  to  2,000  nm 
were  fabricated  on  a  sapphire  substrate  using  electron  beam 
lithography  and  a  lift-off  process  (at  a  fixed  linewidth-to-period 
ratio  of  25%).  An  identical  nickel  nanostructure  was  prepared  on 
a  fused-silica  reference  substrate.  As  discussed  in  Supplementary 
Section  SI,  compared  with  fused  silica  (which  has  an  averaged  mean 
free  path  for  heat-carrying  phonons  of  ~2  nm),  the  averaged  mean 
free  path  of  heat-carrying  phonons  in  sapphire  is  of  the  order  of 
100-150  nm  at  room  temperature.  All  experiments  were  carried  out 
at  room  temperature  and  in  high  vacuum  (~10-6  torr).  The  nickel 
lines  were  impulsively  heated  by  ultrashort  (25  fs)  pulses  from  a 
high-power  Ti:sapphire  laser-amplifier  system.  Pump  pulses  at  a 
wavelength  of  800  nm  were  focused  to  a  fluence  of  2.5mjcm-2, 
in  a  spot  roughly  five  times  the  size  of  the  nanostructured  nickel 
sample  to  ensure  uniform  heating.  Owing  to  the  high  transmission 
of  the  sapphire  and  fused-silica  substrates  at  800  nm,  only  the 
nanostructure  was  heated. 

To  probe  the  heat  flow  dynamics,  an  ultrafast,  coherent  soft 
X-ray  beam  is  diffracted  from  the  nanostructure  at  a  variable 
pump-probe  time  delay.  The  soft  X-ray  beam  is  generated 
through  the  high-order  harmonic  generation  (HHG)  process25  by 
focusing  ~1  mj  pulses  from  the  Ti:sapphire  amplifier  into  a  hollow 


Phonon-transport  regime 


Inverse  Knudsen  number  =  L/A 


Figure  2  |  Analytical  prediction  of  the  Fourier  and  ballistic  components  of 
the  resistivity  for  thermal  transport  away  from  a  half-cylinder  of 
diameter  L.  When  the  linewidth  L  is  much  larger  than  the  phonon  mean 
free  path  A,  the  Fourier  component  (blue)  is  the  largest,  but  for  L<A,  the 
ballistic  resistivity  (red)  is  dominant.  Figure  4  shows  the  first  measurement 
of  the  ballistic  correction  to  the  Fourier  resistivity. 

waveguide  filled  with  argon  gas.  The  resulting  nonlinear  interaction 
between  the  intense  laser  pulses  and  the  atoms  in  the  gas  generates 
odd  harmonics  of  the  driving  laser  frequency.  Absorption  of  the 
argon  gas  and  aluminium  filters  (used  to  block  the  residual  800  nm 
light)  eliminate  all  but  the  25th,  27th  and  29th  harmonic  orders, 
with  a  spectrum  peaked  at  29  nm.  The  full  spatial  coherence  and 
short  wavelength  of  the  soft  X-ray  beams  from  the  waveguide 
make  them  ideal  for  highly  sensitive  interferometric  measurements 
of  nanoscale  surface  deformation26.  Furthermore,  the  very  short 
duration  of  the  HHG  pulses  (sub-10  fs;  ref.  27)  makes  possible 
unprecedented  temporal  resolution  for  dynamics  experiments.  The 
soft  X-ray  harmonics  are  refocused  to  a  ~  1 00  pm  spot  on  the  sample 
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Figure  3  |  Normalized  dynamic  soft  X-ray  diffraction  signal  for  Ni  lines  on 
sapphire,  a-d,  Ni  linewidths  are  810  nm  (a),  350  nm  (b),  190  nm  (c)  and 
80  nm  (d).  The  signal  in  each  case  consists  of  a  sharp  rise  owing  to 
impulsive  laser  heating,  a  thermal  decay  owing  to  interfacial  thermal 
transport  and  an  oscillation  owing  to  surface  acoustic  wave  propagation. 
The  red  and  blue  curves  show  the  fit  using  the  Fourier  law  for  thermal 
transport  in  the  substrate  and  either  large-scale  interfacial  resistance 
(blue)  or  a  modified  effective  interfacial  resistance  ^Effective/  units  of 
10-9  K  m2  W-1),  which  accounts  for  ballistic  effects  from  a  nanoscale 
heat  source  (red).  The  deviation  between  the  two  fits  increases  with 
decreasing  linewidth. 

using  a  grazing-incidence  toroidal  mirror.  The  HHG  beam  reflected 
and  diffracted  from  the  nanostructure  is  then  recorded  using  a  CCD 
(charge-coupled  device)  camera  (Fig.  lb). 

The  change  in  soft  X-ray  probe  diffraction  as  a  function  of  delay 
after  the  infrared  heating  pulse  for  varying  nickel  linewidth  is  shown 
in  Fig.  3.  The  signal  decays  as  the  heat  in  the  nickel  nano -lines 
dissipates  into  the  substrate.  The  rate  of  this  thermal  decay  depends 
on  the  thermal  transport  across  the  nickel/substrate  interface  and  in 
the  substrate.  The  rapid  oscillations  in  the  signal  are  due  to  surface 
acoustic  wave  (SAW)  propagation  in  the  nanostructure.  Periodic 


surface  strain  caused  by  the  patterned  heating  generates  the  SAW; 
we  confirm  that  this  is  the  source  of  the  oscillations  by  comparing 
the  measured  frequency  of  oscillation  with  the  Rayleigh  velocity  for 
SAW  propagation  in  the  substrate28. 

To  interpret  our  results,  we  use  an  analysis  similar  to  that 
used  for  TTR  measurements  of  interfacial  thermal  resistivity.  To 
extract  the  cooling  dynamics  of  the  nickel  nanostructure  from 
data  shown  in  Fig.  3,  we  implemented  a  multiphysics  model 
including  thermal  transport,  thermomechanics  and  Fresnel  optical 
propagation  (see  Supplementary  Section  S2).  The  model  takes 
the  full  two-dimensional  geometry  into  account.  First,  thermal 
transport  across  the  interface  is  governed  by  the  conservation  of 
energy  across  the  surface,  where  x  and  y  are  along  and  away  from 
the  surface  respectively 


PNi^Ni^Ni 


d  TNi(x,f) 
d  t 


(Tm{x,t)-Ts(x,y  =  0,t)) 

^Effective 


where  rEffective  is  the  interface  resistivity  that  will  be  modified  later 
to  include  ballistic  effects  away  from  the  interface,  pNi  and  CNi 
are  the  density  and  volumetric  specific  heat  in  the  nanostructure, 
respectively  and,  TNi  and  Ts  are  the  time-  and  position-dependent 
temperatures  in  the  nanostructure  and  substrate,  respectively.  Heat 
propagation  is  modelled  by  numerically  solving  the  heat  equation 


dTs(x,y,t) 

dt 


G 


V2Ts(x,y,f) 


(1) 


in  the  substrate  beneath  a  single  nickel  line  using  periodic  boundary 
conditions  to  account  for  the  neighbouring  lines.  In  equation  (1), 
ks  and  Cs  are  the  thermal  conductivity  and  volumetric  specific 
heat  of  the  substrate  respectively.  Second,  thermal  expansion  of  the 
substrate  surface  A  hs  is  calculated  from  the  thermoelastic  equation 
assuming  a  stress -free  interface29,  which  involves  a  weighted 
summation  over  the  temperature  distribution  at  each  point 


Ahs(x)  =  2 


(l  +  VsK 

3jt 


If 

J  X\  j  V 


T(xi,y) 


y  dxdy 
(x-Xi  )2+y2 


where  as  is  the  coefficient  of  linear  thermal  expansion  of  the 
substrate  and  vs  is  Poisson’s  ratio.  Thermal  expansion  of  the  nickel 
lines  is  also  accounted  for  using  A/iNi(x)  =  aNirNi(x).  Finally, 
Fresnel  optical  propagation  and  diffraction  from  the  resulting 
surface  profile  gives  the  relative  signal  contributions  from  the 
nickel  lines  and  the  substrate.  We  parameterize  the  results  of  this 
numerical  model  in  terms  of  an  effective  resistivity  rEffective  that  acts 
at  the  interface  and  includes  both  conventional  thermal  boundary 
resistivity  and  size-dependent  ballistic  effects  in  the  substrate  owing 
to  the  confined  heat  source,  and  fit  to  the  time-resolved  thermal 
decay  data.  We  use  the  parameter  resistivity  (r,  units  m2KW_1) 
rather  than  the  resistance  (R,  units  K  W_1),  as  the  latter  depends  on 
the  contact  area  of  the  nanostructured  interface. 

To  qualitatively  understand  the  nature  of  the  ballistic  thermal 
transport,  we  modified  the  model  proposed  by  Chen  for  heat  flow 
from  a  nanoscale  sphere  into  an  infinite  substrate13,  which  sepa¬ 
rately  calculates  the  diffusive  and  ballistic  contributions  to  the  total 
flux17.  We  implemented  the  model  using  a  cylindrical  geometry 
and  extracted  the  resistivity,  which  is  inversely  proportional  to 
the  heat  flux.  Details  of  our  model  are  given  in  Supplementary 
Section  S3.  The  results,  shown  in  Fig.  2,  illustrate  the  nature  of 
the  quasi-ballistic  heat-transfer  effect  as  a  function  of  the  inverse 
Knudsen  number  L/A.  For  linewidths  larger  than  the  mean  free 
path,  the  resistivity  predicted  by  the  Fourier  law  dominates  and  the 
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Figure  4  |  Measured  effective  thermal  resistivity  for  nickel 
nanostructures  of  width  L  deposited  on  fused-silica  and  sapphire 
substrates.  The  blue  and  red  dotted  horizontal  lines  show  the  large-scale 
resistivity  ttbr  for  data  from  the  fused-silica  and  sapphire  substrates.  The 
blue  and  red  dashed  curves  show  model  predictions  for  the  ballistic 
resistivity  correction  r bc  ~  A/L,  with  yips  =  2  nm  and  Asa  =  120  nm, 
respectively.  The  error  bars  indicate  the  standard  deviation  in  the  fits 
to  the  data. 


ballistic  resistivity  is  negligible.  However,  when  the  dimension  of 
the  heated  interface  is  comparable  to  or  smaller  than  the  phonon 
mean  free  path  in  the  substrate,  the  Fourier  law  significantly 
under-predicts  the  resistivity  (or  over-predicts  the  heat  flux).  The 
under-prediction  in  the  Fourier  law  is  given  by  rBaiiistic/ /Fourier  >  so 
we  introduce  a  ballistic  correction  resistivity  rBc  proportional  to 
this  ratio  (and  therefore  proportional  to  A/L),  as  shown  by  the 
red  line  in  Fig.  2.  The  Fourier  resistivity,  rFourie r>  is  automatically 
accounted  for  in  our  numerical  model  for  heat  propagation  in  the 
substrate  (equation  (1)).  We  combine  the  ballistic  correction  with 
the  size-independent  thermal  boundary  resistivity  for  transport 
across  the  interface  to  get  an  effective  resistivity  given  by 


^Effective  —  rTBR  +  ?BC 


In  our  analyses,  we  extract  the  effective  resistivity  /Effective  for  each 
linewidth  from  the  best  fits  (minimized  root  mean  squared  error) 
to  our  data.  By  definition,  large  linewidth  measurements  on  both 
samples  yield  the  large-scale  values  of  thermal  boundary  resistivity 
rTBR  (equivalent  to  TTR  measurements  of  thermal  boundary 
resistance  on  thin  films),  as  shown  by  the  horizontal  dotted  lines 
in  Fig.  4.  The  ballistic  correction  rBC  is  given  by  the  deviation  of  the 
measured  resistivity  at  a  given  linewidth  from  the  resistivity  rTBR 
measured  at  large  linewidths  L.  The  best  fits  to  the  data  shown  in 
Fig.  3  (red  lines)  demonstrate  that  the  measured  interface  resistivity 
increases  as  the  heated  region  gets  smaller.  The  measured  values  of 
/Effective,  Ttbr  and  rBC  for  both  sapphire  and  fused-silica  substrates 
are  shown  in  Fig.  4.  We  measure  a  small  deviation  from  rTBR  for 
the  fused-silica  substrate,  because  the  averaged  phonon  mean  free 
path  in  fused  silica  A¥S  ~  2  nm  is  much  shorter  than  the  smallest 
nickel  linewidths.  However,  in  the  case  of  heat  transfer  from  nanos¬ 
tructures  deposited  on  sapphire  substrates,  the  extracted  resistivity 
at  the  smallest  linewidths  is  more  than  three  times  higher  than  the 
large-scale  value.  The  blue  dotted  and  red  dashed  lines  in  Fig.  4 
show  the  predicted  resistivity  values  /  Effective  (with  rBC  =  /  Effective —  /tbr 
proportional  to  A/L)  based  on  the  analytical  ballistic-diffusive 
model  (Fig.  2).  The  best  fits  to  the  experimental  data  give  estimates 


of  the  average  phonon  mean  free  paths  of  ylFs  =  2  nm  (upper 
limit  of  ~5nm  beyond  which  a  good  fit  is  no  longer  possible) 
and  dSa  =  120  =b  10  nm.  As  discussed  above  and  in  Supplementary 
Section  S4,  these  values  compare  well  to  the  expected  value  of  A. 

Our  measurements  demonstrate  that  the  physical  fundamentals 
of  phonon  transport  from  a  nanoscale  heat  source  can  be 
understood  by  considering  a  minimum  heated  region  model13, 
which  proposes  that  Fourier  transport  can  be  corrected  by  requiring 
the  length  scale  of  the  heat  source  to  be  larger  than  the  phonon  mean 
free  path:  L  ^  L  +  A.  The  effective  resistivity  correction  in  this  case 
is  (1  +  A/L),  which  is  the  experimental  dependence  we  observed, 
shown  as  the  red  dashed  line  in  Fig.  4.  This  implies  that  heat  flow 
from  a  nanoscale  heat  source  smaller  than  the  phonon  mean  free 
path  can  be  thought  of  as  emerging  from  a  larger  spot  corresponding 
to  the  size  of  the  phonon  mean  free  path  in  the  substrate. 

This  is  the  first  direct  and  systematic  measurement  of  quasi- 
ballistic  thermal  transport  effects  from  a  nanoscale  heated  region. 
We  observe  the  breakdown  of  conventional  Fourier  heat  con¬ 
duction  from  interfaces  smaller  than  the  phonon  mean  free  path 
A.  We  find  excellent  agreement  between  the  experimental  data 
and  an  analytical  model  that  determines  the  thermal  transport  by 
summing  the  ballistic  and  diffusive  contributions,  and  demonstrate 
that  the  Fourier  law  can  still  be  used  if  corrected  with  an  extra 
size- dependent  resistance,  proportional  to  the  Knudsen  number 
A/L,  accounting  for  ballistic  effects.  Ballistic  effects  are  expected 
to  get  even  stronger  as  nanostructures  shrink  in  size  or  for  sili¬ 
con  substrates  with  even  longer  mean  free  paths  (ASi  ~  250  nm). 
These  experiments  thus  advance  the  understanding  of  heat- transfer 
fundamentals,  and  are  important  for  the  design  and  manipulation 
of  nanoscale  thermal  transport  in  circuits,  thermoelectrics,  photo- 
voltaics  and  other  structures  of  interest  in  nanotechnology. 
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